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Abstract: In this review, we discuss the zero and finite temperature behavior of 
various bipartite quantum and total correlation measures, the skew information-based 
quantum coherence, and the local quantum uncertainty in the thermal ground state of the 
one-dimensional anisotropic XY model in transverse magnetic field. We compare the ability 
of considered measures to correctly detect or estimate the quantum critical point and the 
non-trivial factorization point possessed by the spin chain. 
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1. Introduction 

In recent years, methods of quantum information science have been extensively used in numerous 
fields of physics, especially in condensed matter theory [1]. One of the most important tools borrowed 
from quantum information theory is based on the concept of entanglement, which is thought to be 
the characteristic trait of quantum mechanics. Entanglement has been investigated in many physical 
settings, and it has been widely considered to be the main resource in most of the quantum information 
processing tasks [2-4]. However, in the last decade, it has been shown that entanglement is not the 
only meaningful correlation present in quantum states, that is separable quantum states can also be 
exploited to provide a quantum mechanical speed-up over classical methods. Quantum discord (QD) 
is the most significant correlation measure that can capture quantum correlations more general than 
entanglement [5-7]. The demonstration of the usefulness of QD as a novel quantum resource has 
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triggered a new line of researeh, which aims to further characterize the quantum correlations beyond 
entanglement in the quantum information community. Following QD, many other quantum and total 
correlation measures have been introduced and widely studied from several perspectives [8-11]. 

Quantum phase transitions (QPT) are sudden changes occurring in ground states of many-body 
quantum systems, which are driven by a control parameter in the Hamiltonian describing the system [12] . 
The key feature of QPTs is that they have their roots in purely quantum fluctuations due to the 
Heisenberg uncertainty relation, in contrast to classical phase transitions, which are only driven by 
thermal fluctuations. In principle, QPTs occur at absolute zero temperature, which is unattainable 
experimentally in practice. However, they are more than a mere theoretical construct, since the signatures 
of QPTs can also be observed at finite and experimentally accessible temperatures. In fact, traces of 
QPTs are detectable when the thermal fluctuations are smaller than the de Broglie wavelength of the 
particles in the system. QPTs are related to the energy level crossings occurring in the ground state of 
the many-body system, which result in the discontinuities in the ground state energy. The order of the 
phase transition is then determined based on this discontinuity. While a discontinuity in the ground state 
energy signals a first-order transition, a discontinuity in the first derivative of the ground state energy can 
be observed when the second-order phase transition takes place. 

Since quantum spin chains possess several different kinds of QPTs, they are natural candidates for 
studying quantum critical phenomena in many-body systems. On the other hand, another peculiar 
property of spin chains in a transverse magnetic field is the factorization phenomena [13]. Although 
the ground states of such systems are in general entangled, for some specific Hamiltonian parameters, 
the ground state becomes completely factorized. This phenomenon has its roots in the symmetry of the 
ground state of the system under consideration. The detection and understanding of this phenomena are 
rather non-trivial and have been the subject of many works [14-22]. 

In this review, we intend to cover the literature on the behavior of correlations in various quantum spin 
chains, both near QPT and the factorization point (FP). We summarize the results on the key contributions 
in the field in detail and discuss their implications both at zero [23-59] and finite temperatures [60-70]. 
The remainder of the review focuses on some specific results concerning the correlations, coherence 
and uncertainty in the ID anisotropic XY model in the transverse field. In particular, we summarize 
the analysis of the behavior of four quantum and total correlation measures at zero and finite 
temperatures [65,68]. In order to quantify total correlations, a general measure of non-locality [9] and a 
measure based on Wigner-Yanase skew information [10] are used. Quantum correlations are quantified 
by an experimentally lower bound of the geometric version of QD [11], and entanglement is quantified 
by concurrence [3,4]. On the other hand, coherence and uncertainty are quantified by the regular and an 
optimized version of Wigner-Yanase skew information, respectively [7 1 ,72] . 

This review is organized as follows. In Section 2, the spin-1/2 anisotropic XY model in a transverse 
magnetic field is introduced. QPTs, factorization phenomena and the symmetries of the ground state 
are discussed. Section 3 introduces the quantum and total correlation measures used in this review 
together with their physical interpretations. Section 4 begins with a comprehensive review of the works 
on quantum correlations in quantum spin chains and continues with some detailed specific results on the 
XY model in a transverse magnetic field. Section 5 includes our conclusion. 
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2. Spin-1/2 Anisotropic XY Chain in a Transverse Field 


The Hamiltonian of the one-dimensional anisotropic spin-1/2 XY chain in a transverse magnetic field 
is given by: 


, N N 

^ + (1 - 

i=i i=i 

where are the usual Pauli operators at the j-th site, A denotes the strength of the inverse magnetic 
field, 7 G [0,1] is the anisotropy parameter and N is the number of spins. The Hamiltonian, which is 
in the Ising universality class for 7 > 0, reduces to the XX chain for 7 = 0 and to the Ising model 
in transverse field for 7 = 1. The system possesses a second order QPT in the parameter interval 
0 < 7 < 1 at the critical external field Ac = 1 that separates a paramagnetic (disordered) phase from 
a ferromagnetic (ordered) phase. There is another QPT in the region A > 1 at 7 c = 0, which is of the 
Berezinskii-Kosterlitz-Thouless type and separates a ferromagnet ordered phase in the x-direction from 
a ferromagnet ordered phase in the y-direction. However, this transition will not be deeply explored in 
this review. 

Although the ground state of the considered model is entangled in general, there exists a non-trivial 
factorization line corresponding to 7 ^ -f A“^ = 1. Therefore, the ground state becomes completely 
factorized at any field point satisfying the equation. 


- 


1 

Vl - 7 ^ 


( 1 ) 


The diagonalization procedure of the XY model includes the well-established techniques of 
Jordan-Wigner and Bogoliubov transformations [73,74]. The model Hamiltonian is invariant under a 
TT rotation around the z-axis of all spins (a parity transformation), which can also be stated by saying 
that the system exhibits global Z 2 symmetry. Taking this symmetry into account together with the 
translational invariance of the system, the reduced density matrix of two spins have an X-shaped form, 
and in terms of the magnetization and two-spin correlation functions, it can be written as: 
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Here, r denotes the distance between the spins, and the two-spin reduced density matrix is only 
dependent on the distance between the spins, r = \i — j\, due to the translational invariance. The 
magnetization and two-spin correlation functions are given as [73,74]: 
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where the function G^ is given as follows: 
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and = a/( 7 A sin 0 )^ + (1 + A cos 0)^/2 with (3 = 1 /kT being the inverse temperature. Tracing out 
one of the spins, we can write the single-spin reduced density matrix as follows: 


Po — Pi 


1(1 + (a^) 0 A 


(3) 


where (a^) is the transverse magnetization, and the density matrix is written on the basis of the 
eigenvectors of a^. Note that due to the translational invariance of our system, it is not important which 
spin we trace out. Indeed, all of the single-spin reduced density matrices are the same in this case. 

The paramagnet-ferromagnet transition at Ac = 1 in the 0 < 7 < 1 region has important consequences 
regarding the ground state of the system. At the CPof the QPT, Z 2 symmetry is broken with the 
expectation value of the magnetization in the x-direction, (Sx), being the order parameter for this 
transition with a non-zero value for A > 1, and the ground state of the system becomes two-fold 
degenerate. In real physical settings, the system chooses one of the degenerate ground states as the 
real ground state due to some small external perturbation, which is called the mechanism of spontaneous 
symmetry breaking (SSB). However, in this case, the ground state of the system does not possess the 
symmetries of the Hamiltonian, for example the Z 2 symmetry in our case, and the reduced density 
matrix given in Equation (2) no longer has the X-shaped form. In the vast majority of the studies 
present in the literature, the effect of the broken symmetry is not taken into account, apart from a few 
exceptions [69,70,75-78], due to the complications it introduces into the calculations. We should note 
that, throughout this paper, we also neglect the effects of SSB and consider the so-called thermal ground 
state of the system. The thermal ground state is an equal mixture of these two degenerate states. Indeed, 
it is nothing but the limit fd ^ 00 of the canonical ensemble. 
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where Z is the partition funetion. It is important to note that the thermal ground state is aetually the same 
as the real ground state when there is no ground state degeneraey. Therefore, in the disordered phase, we 
are working with the real ground state of the system. 

3. Correlations, Coherence and Uncertainty 

In this section, we introduce the measures of correlations, coherence and uncertainty, which are 
central to this work. In the following sections, we will extensively investigate the behavior of the 
measures introduced here in the context of criticality and factorization in the anisotropic XY model. 

3.1. Geometric Measure of Quantum Discord 

Quantum discord is the first and the most popular quantum correlation measure in the zoo of measures 
that are more general than entanglement [6,7]. It has been defined in terms of the discrepancy between 
the two classically equivalent generalizations of the quantum mutual information and can be present in 
states that contain no entanglement. 

The amount of total correlations in a given quantum system can be quantified by the mutual 
information: 

I{p^^) = S{p^) + S{p^)-S{p^% (5) 

where p^^ is the density matrix of the total system, and p^ are the reduced density matrices of the 
subsystems and S{p) = —trplog 2 P is the von Neumann entropy. On the other hand, it is possible to 
quantify the classical correlations contained in a quantum system as follows [7]: 

Cip'^) = s(p“) - minY^ptS(pl), (6) 

{n|} V 

where ( 11 ^) is the set of most general quantum measurements performed on subsystem b and p% = 
(/ (g) (g) Y{\)/pk are the post-measurement states of subsystem a after obtaining the outcome k 

with probability pk = tr(J“ 0 from the measurements made on subsystem b. Quantum discord 

is defined as the difference between the total and classical correlations, which, as a result, quantifies the 
quantum correlations: 

6{p^’’)=I{p^’’)-C{p^^). (7) 

However, calculation of quantum discord is a rather difficult task due to the complex optimization 
procedure in the evaluation of the classical correlations, C. In order to overcome this difficulty, a 
geometric measure of quantum discord (GMQD) has been introduced [ 8 ], which measures the distance 
between the given state and the set of zero-discord states. Mathematically, it can be expressed as follows: 

DG(p“')=2min||p“'-xir, (8) 

X 

where the minimum is taken over the set of zero-discord states, and ||.||^ denotes the square of the 
Hilbert-Schmidt norm. 
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The representation of a general bipartite state in a Bloeh basis ean always be written as: 
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where the matriees {Xi : i = 0,1, • • • , — 1} and {1^ : j = 0,1, • • • , — 1}, satisfying tr(XfcX;) = 

tr(yA:^z) = Ski, define an orthonormal Hermitian operator basis assoeiated with the subsystems a and b, 
respeetively. The eomponents of the loeal Bloeh veetors x = {xi}, y = {yj} and the eorrelation matrix 
T = tij ean be obtained as: 


X, = trp“'(X,0/')/v^, 

Vj = ® Yj)/Vm, 

Uj = trp'^\Xi^Yj). ( 10 ) 


Following the definitions made above, reeently, an interesting analytieal formula for the GMQD of an 
arbitrary two-qubit state has been introdueed [11]: 


= 2(trS' - max{ki}), 


( 11 ) 


where S = xx'' + TT* and 
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( 12 ) 


with {ai} = {0,27r,47r} and 6 = arccos{(2trS'^ — QtrS'trS'^ + 9trS'^)-^2/(3trS'2 — (trS')^)^}. 
Furthermore, observing that cos reaches its maximum for Oj = 0 and choosing 9 to be zero, 

a very tight lower bound to the GMQD can be obtained, and it is given by: 

Q(p“'’) = ^(2tr^ - V6tr^2 _ 2(tr^)2). (13) 

o 

We will refer to this quantity as the observable measure of quantum correlations (OMQC). It satisfies 
all of the criteria to be a meaningful measure of quantum correlations on its own. Moreover, it has the 
important advantage of requiring no optimization in the evaluation process. On top of being easier to 
calculate than the GMQD, it can be measured by performing seven local projections on up to four copies 
of the state. Thus, Q{p) is also very experimentally friendly, since one does not need to perform a full 
tomography of the state. 


3.2. Measurement Induced Nonlocality 

After the introduction of the quantum discord, various quantum and total correlation measures 
have been proposed to quantify different kinds of correlations contained in a quantum system. 
Measurement-induced nonlocality (MIN) is a correlation measure that encapsulates the nonlocality 
contained in a bipartite system, in a more general way than the violation of Bell inequalities or 
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entanglement, from a geometrie perspeetive based on von Neumann measurements [9]. The physieal 
setting behind MIN is as follows: We perform von Neumann measurements on one of the subsystems 
of a bipartite quantum system and interpret the differenee of the pre- and post-measurement states as a 
quantifier of nonloeality. The measurement operators are ehosen so that they do not disturb the applied 
subsystem, in order to quantify solely the nonloeal effeet of the measurement. Mathematieally, it is 
defined by (taking into aeeount the normalization): 

= 2max||p“'’ -n“(p“'’)||2, (14) 

where the maximum is taken over the von Neumann measurements 11“ = {H^} that do not ehange p“ 
loeally, meaning n^p“n^ = p“, and ||.||^ denotes the square of the Hilbert-Sehmidt norm. Although 
MIN, as given by Equation (14), has no elosed form formula for an arbitrary bipartite state, it ean be 
ealeulated analytieally for pure states of arbitrary dimension and for 2 x N dimensional mixed states. 
MIN for a 2 X 2 dimensional system (two-spin system), whieh will be the foeus of this work, ean be 
analytieally evaluated as: 

f 2(trTT* — if x 0, 

N(p)={ IWP ^ (15) 

l2(trrr*-A3) iff=0, 

where TT* is a 3 x 3 dimensional matrix with A 3 being its minimum eigenvalue, and ||a; IP = 
with X = {xi,X 2 , x^y. Due to the speeial form of the density matrix eonsidered in this work, whieh is 
given in Equation (2), and sinee the loeal Bloeh veetor x is never zero in our investigation, MIN takes 
the simple form: 

^(P) =4(pL + pL)- (16) 


3.3. Wigner-Yanase Skew Information 

Wigner-Yanase skew information (WYSI) mainly quantifies the information eneapsulated in a 
quantum state with respeet to an observable that is not eommuting with it [79] . In reeent years, WYSI has 
attraeted a lot of attention in very different eontexts, leading to different physieal interpretations [80-83]. 
The mathematieal expression for WYSI is given as: 

/(p,/q = -fTr[VP.Af, (17) 

where K denotes an observable, i.e., a self-adjoint matrix. A very important property of WYSI is that 
it eaptures the genuine quantum uneertainty of a given observable in a eertain quantum state. In other 
words, it does not get affeeted by the elassieal mixing. The basie properties of the WYSI ean be listed 
as follows: 

• It is upper bounded by the varianee with respeet to the eonsidered observable, I{p,K) < 
V (p, K) = TrpiT^ — (TrpiT)^, with the equality reaehed for pure states. 

• /(p, K) is eonvex, sueh that: 

/(5^A,p„iT) <5 ^A,J(p,,JO, 


( 18 ) 
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where pi is an arbitrary quantum state and A* is a probability distribution satisfying the eonstraint 

Zli = 1- 

• For a given bipartite state p°^ and an observable K aeting on the subsystem a, one has: 

l{p^\K®l)>l{p\K). (19) 

Coherenee is one of the most signifieant properties that separates quantum states from elassieal 
ones. Although various different intuitive methods to measure the coherenee of a given quantum state 
are known, there were no general frameworks for quantifying coherence. Recently, a list of defining 
properties for the coherence measures has been introduced in order to properly quantify coherence [84] . 
These properties can be stated as follows: a coherence monotone is: (i) zero if and only if the state 
is incoherent; (ii) non-increasing under incoherent operations; and (iii) non-increasing under classical 
mixing of quantum states. It has been shown that WYSI satisfies all of the criteria for coherence 
monotones [71] and is defined as a JT-coherence measure, meaning that the quantifier of coherence 
is contained in p when measuring the observable K. In fact, the intuitive physical explanation behind 
this can be expressed as follows: WYSI gives information only about the quantum uncertainty appearing 
as a result of a measurement, and this uncertainty has its root directly in the coherence embodied in a 
state. Here, we follow this line of work and adopt WYSI as a coherence measure, referring to it as local 
quantum coherence (LQC). 

Furthermore, it is possible to obtain a lower bound for WYSI by dropping the square root on the 
density matrix: 

Ipp,K) = -jJ’c[p,K]\ ( 20 ) 

This lower bound is particularly important due to its easy experimental accessibility. It can be 
measured using two programmable measurements on an ancillary qubit. Note that LQC can be calculated 
both for single and bipartite systems. Naturally, for bipartite systems, it is written as I{pab, Ka ® Ib), 
which quantifies the coherence with respect to a local observable acting on the first subsystem. 

We have already mentioned that WYSI does not get affected by the classical mixedness of the 
quantum state. Motivated by this property of WYSI, in [72], the local quantum uncertainty (LQU) 
of a given state with respect to an observable is defined as the minimum achievable WYSI on a single 
local measurement: 

U^^ = mmI{p,K^^), (21) 

where F denotes the spectrum of K^, and the minimization over a chosen spectrum of observables leads 
to a specific measure from the family. However, for a two-qubit system, all of the members of the family 
turn out to be equivalent. Then, the LQU can be analytically calculated as: 

Ua{PAb) = 1 — AmaxIkkAs}, 

where A^ax is the maximum eigenvalue of the 3x3 symmetric matrix Wab^ whose elements are given by: 

{WAB)ij = Tr {^JpAB{criA ® Ib)\/ PAB{crjA ® Ib)} , 




Entropy 2014 , xx 


9 


where indiees i,j = {x,y,z} are given for the usual Pauli operators. We note that Equation (21) is 
normalized to one for maximally entangled pure states and, moreover, reduces to the linear entropy for 
any pure bipartite state. 

Lastly, we introduce a measure for total correlations based on the WYSI, which we will refer to as the 
Wigner-Yanase skew information-based measure (WYSIM) [10]. As I{p, X) depends both on the state 
p and the observable X, an averaged quantity has been introduced in order to get an intrinsic expression: 

Q(p) = 5^/(p,X,), (22) 

i 

where {Xi} is a family of observables, which constitutes an orthonormal basis. It is possible to use 
Q{p) to capture the total information content of a bipartite quantum system with respect to the local 
observables of the subsystem a as follows: 

(23) 

i 

which is independent of the choice of a specific orthonormal basis {Xi}. As a result, the difference 
between the information content of the composite system and the Kronecker product of the local 
subsystems 0 p^ with respect to the local observables of the subsystem a can be introduced as a 
measure for total correlations for 

= 5(Q.(p“‘) - Q«(p“)). (24) 

where we add a normalization factor 2/3. 

Although many entanglement and other correlation measures mostly involve complicated 
optimization procedures, WYSI enjoys the advantage of being able to be calculated quite 
straightforwardly. Another entropic measure for total correlations, which is also easy to calculate, is 
the quantum mutual information introduced in Equation (5), which has also been shown to detect the 
CPs of various quantum critical systems. At this point, it is important to note that quantum mutual 
information and WYSI have different physical interpretations and reveal different information about 
the spin chain under consideration, which, in return, extends our understanding about the behavior of 
quantum information in such systems. 

3.4. Concurrence 

In order to quantify the entanglement content of a two-qubit density matrix, we utilize concurrence, 
which is the most common way to measure entanglement for these systems [3,4]. Eirst, we need to 
calculate the time-reversed or spin-flipped density matrix p, which is given by: 

p = ((T^ (8) a^)p*{a'^ ® a^). 

Here, a'^ is the Pauli spin operator and p* is obtained from p via complex conjugation. Then, concurrence 
reads: 


C{p) = max |o, \/~Xi - - \/^ - |, 


(25) 
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where {Aj} are the eigenvalues of the produet matrix pp in decreasing order. For the simple form of the 
reduced density matrix given in Equation (2), concurrence reduces to: 

C = 2max{0, jpu] - \p 22 l \p 2 z\ - yjpiipii}- (26) 


4. Behavior of Correlations 

The investigation of various quantum correlation measures in the ground states of spin models have 
been an active area of research in the last decade. Before this line of research was initiated, Preskill 
argued that better understanding of entanglement in strongly-interacting many-body systems may deepen 
our knowledge of these systems [85]. Following this idea, one of the first steps was taken in the direction 
of analyzing the single-site and two-site entanglement in the anisotropic XY chain in the transverse 
magnetic field (and in a limit of the same model, the transverse Ising chain) [23]. In this work, along 
with many important results, it has been shown that the quantum critical point in these systems can 
be identified by looking at the behavior of the entanglement in the ground state of the model, as 
quantified by concurrence and von Neumann entropy. Both measures are calculated for nearest-neighbor 
or next nearest-neighbor spins, and it has been concluded that the bipartite entanglement close to the 
CP in these systems is short-ranged. Moreover, the study includes the exploration of the bipartite 
entanglement in finite temperatures. It can also be seen that the entanglement is not very susceptible 
to the effect of temperature. In [24], entanglement in the finite-sized XY chain is investigated. It has 
been shown that it obeys a scaling behavior near the critical point of the QPT in the model. In order 
to understand the connection between the non-analyticities of the bipartite entanglement measures and 
the QPTs, a general framework has been introduced. In [25], it has been discussed that the behavior of 
bipartite entanglement measures is directly related to discontinuities occurring in the ground state energy. 
Therefore, the conditions to detect QPTs are obtained under the assumptions that these discontinuities are 
not accidentally canceled or no artificial discontinuity is created in the optimization procedure involved 
in the calculation of the measure. The last assumption is particularly important, because there are some 
explicit examples in the literature where such cases occur [34]. 

On the other hand, the behavior of multipartite entanglement in spin chains has also been investigated. 
In [37-39], based on a multipartite entanglement measure, which is a generalization of a global 
entanglement measure [86], it has been demonstrated that it is also possible to find the signatures of the 
QPT in some quantum critical models. Moreover, the multipartite entanglement has been determined 
to be larger than the amount of bipartite entanglement in the vicinity of QPT, suggesting that the 
entanglement is distributed over the whole lattice near criticality. The authors also have extended 
the framework introduced in [25] to the case of multipartite entanglement. Recently, in [40] and in 
[41], multiparticle entanglement in anisotropic XY model in transverse field have been explored using 
different criteria for detecting the entanglement. Both papers show that in spite of the lack of bipartite 
entanglement, multipartite entanglement is present nearly for all values of the external field and the 
anisotropy parameter. Furthermore, they show that both entanglement criteria obey the finite-size scaling 
behavior of the XY model. 

We have already mentioned that the symmetry breaking in the ground state of the considered spin 
chain models is generally not taken into account in the majority of the studies on this subject. However, 
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there are several works exploring the effects of SSB on entanglement near criticality. The initial work 
in this direction was made in [75], where it has been shown that the symmetry breaking does not 
effect entanglement in the ground state of the XXZand Ising model in general. Indeed, through an 
inequality involving the correlation functions of the system introduced in this work, it has been possible 
to determine whether the symmetry breaking will change the amount of entanglement. Later, in [76], the 
work of [75] was extended to the case of the XY model in transverse field and the Lipkin-Meshkov-Glick 
model, and it has been demonstrated that for both of these models, there has been an enhancement in the 
entanglement when the symmetry breaking effect has been taken into account. Furthermore, multipartite 
entanglement has been shown to be affected more severely than the bipartite entanglement by the SSB 
near the quantum critical point [77,78]. 

Many important studies followed this line of work, investigating the interplay between entanglement 
and QPTs in different spin chain models, such as XXZ, XYZ, XY with the Dzyaloshinskii-Moriya 
(DM) interaction, Lipkin-Meshkov-Glick (LMG), etc., both in the thermodynamic limit and in the 
few-body cases. 

Apart from the relation between the QPTs and entanglement, another important question is the 
dynamics of entanglement in the quantum spin chain, which was also addressed [29]. Here, the authors 
analyze how an initially entangled bipartite state evolves under the dynamics governed by the XY 
Hamiltonian as a function of the system parameters and the distance between the initially prepared spins. 

Now, we turn our attention to the relationship between QPTs and quantum correlation measures that 
are more general than entanglement. The interest in this subject has increased after the introduction of 
quantum discord, which is the first and most widely explored measure of quantum correlations beyond 
entanglement. Following quantum discord, many other measures have been introduced, some of which 
we will also explore in the context of the XY model. 

The pioneering work on the investigation of general quantum correlations in spin chains was made 
in [46] . In this work, the author considered the transverse Ising chain together with the antiferromagnetic 
XXZ chain and analyzed the behavior of quantum discord for two nearest-neighbor spins, close to the CP 
of both systems. It has been shown that the quantum discord gets larger close to the CP for both of the 
considered models, but does not become maximum in the Ising model case. Due to this abrupt increase, 
the derivative of quantum discord is able to signal the presence of the QPT. Following this, a similar 
discussion to the one made in [25] has also been considered, and the relation between the discontinuities 
in the reduced density matrix elements and the quantum discord has been put forward. On the other 
hand, while the classical correlations decrease as one gets close to the QPT for the XXZ chain, they 
show a monotonically increasing behavior in the Ising model. A generalization of this work has been 
made in [47], which gives a general recipe for the calculation of quantum discord in states possessing 
Z 2 symmetry, a general spin flip symmetry exhibited by most of the spin chain Hamiltonians, but not 
necessarily exhibited by the ground state density matrix in the ordered phase, as we have mentioned. In 
order to illustrate the findings, the author has investigated the behavior of correlations in XXZ, Ising and 
LMG models, both in the thermodynamic limit and in the finite size cases. He has shown that the effects 
of first-, second- and infinite-order QPTs in these systems can be found in the classical correlations 
and quantum discord, either via the measure’s or its derivatives behaviors near criticality. Furthermore, 
the scaling law obeyed by the derivative of the quantum discord was shown to differ compared to that 
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of the entanglement. A eomprehensive study of the two-site sealing of quantum diseord in the XXZ, 
XY and Ising models in a transverse magnetie field has been done in [59]. It is important to note that 
the geometrie version of the quantum diseord has also been explored in the anisotropie XY model [57]. 
It has been shown that the derivative of the geometrie quantum diseord is also singular at the CP, together 
showing a universal finite-sized sealing behavior. 

The analysis of quantum eorrelations at finite temperatures has also attraeted mueh attention due to 
the eontrast of the obtained results with respeet to the ones obtained for entanglement measures. One 
of the major eontributions on this subjeet has been made in [60-62], in whieh the quantum diseord in 
the anisotropie XY, XXZ and Ising models in a transverse magnetie field has been explored at both 
zero and finite temperatures. They have shown that through the behavior of bipartite, nearest-neighbor 
thermal quantum diseord (TQD), whieh is merely the quantum diseord ealeulated at finite temperature, 
it is possible to deteet the position of the QPT in the parameter spaee of the model Hamiltonian at 
finite temperatures. The partieular importanee of this work also stems from the faet that entanglement 
was not able to signal the CP at non-zero temperatures, as well as the thermodynamie quantities, sueh 
as entropy, speeifie heat, magnetization and magnetie suseeptibility. Sinee the exaet zero temperature 
is not attainable in real physieal systems, sueh a result proves the advantage of quantum diseord in 
experimental systems. Furthermore, the ability in estimating the eorreet values of CPs (the ones at 
zero temperature) in these systems was eompared using the quantum diseord, entanglement and other 
thermodynamie quantities. It has been revealed that although at temperatures very elose to zero, all of 
these quantities pinpoint the CPs with aeeeptable sueeess, as the temperature inereases, quantum diseord 
outperforms other measures. In another work [50], the TQD has also been investigated in two limits 
of the anisotropie XY model in a transverse field; the XX and Ising models. The authors have shown 
that the quantum eorrelations in the system ean inerease with temperature in high magnetie fields. They 
suggest that the explanation behind this effeet is: when the magnetie field is on, the low-lying exeited 
states are more eorrelated than the ground state, and with inereasing temperature, the system favors those 
exeited states. 

The investigation of eorrelations for distant bipartite spin parrs is also an important subtopie in the 
investigation of eorrelations in spin ehains. We have already mentioned that the bipartite entanglement 
near eritieality in the spin models is generally short ranged; beyond next nearest neighbors, it is not 
possible to find a highly entangled pair. However, this is not the ease for quantum eorrelations that are 
more general than entanglement. In [50,54], it has been shown that in the XY model, in eontrast to the 
short-ranged entanglement, as the system approaehes the ordered phase, a long-ranged quantum diseord 
develops in the system. A more systematie approaeh to the long-range eorrelations in the same model 
has been exhibited in [67], where both zero and finite temperature eases have been explored. For the zero 
temperature ease, it has been demonstrated that the quantum diseord of two-spins, whieh are 15 spins 
apart, was able to eapture the QPT. Indeed, it behaves like an order parameter of the transition, being 
zero in the disordered phase, while beeoming non-zero with a finite jump in the ordered phase. On the 
other hand, at finite temperatures, quantum diseord still seems to be able to spotlight the QPT up to 15 
spin separations; however, it beeomes very fragile to the effeets of temperature. 

Similar to the ease of multipartite entanglement, there has also been an interest in the definition, 
eharaeterization and operational meaning of the QD generalized to multipartite states, whieh goes by 
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the name, global quantum diseord (GQD). A widely-accepted definition of GQD has been made in [52] 
by a systematic extension of the original definition of bipartite QD, which is done by rearranging it in 
terms of relative entropy and local von Neumann measurements. The authors have also demonstrated 
the behavior and content of GQD in Werner-GHZstates and the Ashkin-Teller spin chain. The latter 
possesses an infinite-order QPT in its phase diagram. It has been shown that GQD is able to capture the 
presence of this CP where bipartite QD is not able to do so. The results on GQD are extended to various 
spin chain models at both zero and finite temperatures [64,65]. The behavior of GQD has been studied 
in finite-sized transverse Ising, cluster-Ising and XX models in a transverse magnetic field with open 
boundary conditions at finite temperature, and it was shown that the non-trivial changes in these systems 
can be captured by looking at the GQD. Furthermore, it has also been demonstrated that, for the Ising 
model, finite-sized scaling of the GQD is characterized by the universal critical exponents belonging to 
the Ising universality class. 

The factorization phenomena that we introduce in Section 2 is not always so straightforward to detect 
with the help of the correlation measures. First of all, the point where the ground state of the system 
factorizes is always in the ferromagnetic (ordered) phase of the system, as can be seen from Equation (1). 
Thus, in principle, in order to see a clear signature of the factorization on the correlation measures, 
one needs to consider the effects of symmetry breaking in the ground state of the system. However, 
even when SSB is ignored, bipartite entanglement measures (concurrence) for a fixed distance between 
two spins are able to detect it by being exactly equal to zero at this point. On the other hand, QD 
for the thermal ground state of the XY model can show the existence of this factorized point by the 
single crossing of the QD curves for different bipartite spin distances, i.e., by having the same value 
independent of the spin-spin distance [69]. In a recent work [67], a different approach has been taken in 
order to understand the nature of the factorization phenomena. The authors considered a finite-sized spin 
chain, looked at the energy spectrum and saw that right at the factorization field, A/, the energy levels of 
the ground and the first excited states cross, which are of opposite parity. What happens at this point is 
actually a transition between the different parity ground states. They have also examined the factorization 
phenomena at finite temperatures and seen that as the temperature increases, the factorization point 
widens and becomes a region of separability. 

We have given an overview of the main results on the behavior of QD in different spin chains. 
However, there are many other important works done on this topic, such as examining the nonlocality 
and other measures of non-classicality in various spin models with more exotic interactions in both 
few-body systems and in the thermodynamic limit. An emerging subject in this field is to examine the 
correlations in systems possessing a topological phase transition. Topological phase transitions cannot 
be characterized by an order parameter or a broken symmetry in the system, therefore laying outside the 
scope of the traditional way of understanding phase transitions [87]. Thus, they are characterized by the 
change in some topological properties in the ground state of the system. Attacking this problem with 
the tools of quantum information theory has also given some important insights about the nature of the 
topological phase transitions [88-93]. 

In the following sections, we will extensively review our contributions to the investigation of 
correlation measures in spin chains, which involves the study of the correlation, coherence and 
uncertainty measures introduced in Section 3 in the one-dimensional anisotropic XY model in a 
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transverse magnetie field, both at zero and finite temperatures, with a speeial emphasis on the 
deteetion of the faetorization phenomena and ambiguities eaused by the optimization procedure of the 
correlation measures. 


4.1. Correlation Measures 


We begin with the discussion of the total correlation measures as quantified by WYSIM and MIN, 
which are, together with their derivatives, presented in Figure 1. The measures are calculated for nearest 
neighbors and plotted as a function of the external field A for kT = 0, 0.1, 0.5 and 7 = 0.001, 0.5, 1. It 
can be seen that both measures behave in a very similar fashion for 7 = 1 and 7 = 0.5 at all temperatures. 
Their behavior differs in the case of 7 = 0.001, where WYSIM makes a sharper jump to a finite value 
than MIN right at the CP. We can also see that the correlation content of the two-spin reduced density 
matrix increases as a function of A, meaning that the system is more correlated in the ferromagnetic 
phase. Right at the CP Ac = 1, both measures make a finite jump, resulting in a divergence in their first 
derivatives, which points to the existence of the second-order QPT of the considered system. 


Figure 1. Left panel The thermal total correlations quantified by measurement-induced 
nonlocality (MIN) and the Wigner-Yanase skew information-based measure (WYSIM) as 
a function of A for 7 = 0.001, 0.5, 1 at kT = 0 (solid line), kT = 0.1 (dashed line) and 
kT = 0.5 (dotted line). Right panel The first derivatives of MIN and WYSIM as a function 
of A for 7 = 0.001, 0.5, 1 at kT = 0 (solid line), kT = 0.1 (dashed line) and kT = 0.5 
(dotted line). The graphs are for the first nearest neighbors. 
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We can now turn our attention to the investigation of quantum correlations as quantified by OMQC 
and concurrence, and their derivatives, presented in Figure 2. As in the total correlation case, the graphs 
are for nearest neighbors and plotted as a function of the external field A for kT = 0, 0.1, 0.5 and 




































Entropy 2014 , xx 


15 


7 = 0.001, 0.5, 1. Again, we can see that both measures are capable of detecting the CP in the system 
though the divergence in the first derivative. It can be seen that concurrence is not very susceptible to 
the effects of temperature, since it remains zero nearly for all parameter values at high temperatures. 
We also note for 7 = 0.001 that the behavior of OMQC and concurrence is very similar to that of MIN 
and WYSIM, respectively. This suggests that at this parameter range, the correlations contained in the 
system are mainly quantum. 


Figure 2. Left panel The thermal quantum correlations quantified by observable measure of 
quantum correlations (OMQC) and concurrence as a function of A for 7 = 0.001, 0.5, 1 at 
fcT = 0 (solid line), kT = 0.1 (dashed line) and kT = 0.5 (dotted line). Right panel The 
first derivatives of OMQC and concurrence as a function of A for 7 =0.001, 0.5, 1 at kT = 0 
(solid line), kT = 0.1 (dashed line) and kT = 0.5 (dotted line). The graphs are for first 
nearest neighbors. 



(a) 


(b) 


Although the divergence of the derivatives of the measures is clearly pronounced for all values of 
7 at kT = 0, as the temperature increases, the finite jump of the measures near QPT smooths out, 
causing the divergences in the derivative to disappear. This is an expected result, since, as we have 
mentioned earlier, as temperature increases, thermal effects dominate the genuine quantum fluctuations. 
Additionally, similar results can be obtained if we look at the correlations for next nearest neighbors; 
however, they provide no additional information about the system, so we do not present them here. 

Focusing on the signatures of the factorization phenomena, we first observe that for 7 = 0.5, 
concurrence becomes exactly zero at the factorization field A/ ~ 1.15, which can be calculated from 
Equation (1). More interestingly, in all of the remaining three correlation measures considered in this 
work, only WYSIM is able to detect the existence of the factorized ground state by a non-analytic 
behavior in its derivative. This is a rather peculiar property of WYSIM, since, as mentioned, we do 
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not take into account the effects of SSB. No other correlation measure in the literature is capable of 
revealing this phenomena without explicitly considering the effects of SSB with a single evaluation of 
the measure. The intersection point of the QD calculated using the thermal (symmetry protected) state for 
different spin distances can detect factorization in this model; however, for WYSIM, a single calculation 
is sufficient. 

In Figure 3, we present our results on the performance of the measures in correctly estimating the CP 
at finite temperatures. We know that the divergent first derivative disappears and becomes a smoother 
function as the temperature increases. However, at sufficiently low temperatures, it is still possible to 
extract some information from the behaviors of the correlation measures in the vicinity of the QPT. The 
strategy we follow to estimate the CP at a finite temperature is to look for a local maximum or minimum 
at the first derivative and identify this point as the estimated CP. This extremum point can be identified 
easily by checking the second derivative of the measure, as is done here. A first glance. Figure 3 suggests 
that the accurateness of the estimation is strongly dependent on the parameters of the Hamiltonian. 


Figure 3. The estimated values of the CPas a function of kT for three different values 
of the anisotropy parameter 7 = 0.001, 0.5, 1. The CPs in the graphs are estimated by 
OMQC (denoted by o), WYSIM (denoted by +), MIN (denoted by *) and concurrence 
(denoted by x). Concurrence is not included for 7 = 1 and r = 2, since it vanishes at 
even very low temperatures. 
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The graphs in the upper panel of Figure 3 are plotted for the nearest-neighbor spins for 7 = 0.001, 
0.5, 1 as a function of temperature. The considered temperature range is from kT = 0 to kT = 0.1. 
In the case of 7 = 1, we see that the best and the worst estimators turn out to be the concurrence 
and WYSIM, respectively, with low relative error between them. OMQC and MIN nearly estimate 
the same value for CP, with MIN outperforming OMQC only very slightly as we approach kT = 0.1. 
For 7 = 0.5, as OMQC shows a remarkable accuracy by approximately predicting the exact value of the 
CP at kT = 0, MIN deviates from the true value rapidly and by a large amount. Lastly, when we change 
the anisotropy parameter to be 7 = 0.001, we observe that OMQC and MIN behave in the exact same 
manner and perform very poorly in the estimation process, while both concurrence and WYSIM predict 
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the position of the CP better, with WYSIM being the closest to the original value. On the other hand, in 
the lower panel of Figure 3, the next nearest-neighbor spins for the same values of 7 and the same range 
of temperature as the upper panel are plotted. Again, starting with the the 7 = 1 case, we notice that 
WYSIM, MIN and OMQC lose their accuracy for the same amount, and the concurrence is not plotted 
for this case, since it quickly decays to zero, even for very small values of kT. Note that in all cases 
considered in Figure 3, only for these values of parameters do we see a measure, OMQC, estimating 
the CP as less than its true value. Considering the 7 = 0.5, the total correlation measures (WYSIM 
and MIN) perform worse than the quantum correlation measures (OMQC concurrence). Finally, for 
7 = 0.001, we observe a very similar pattern with the graph plotted for the nearest neighbors. 

A similar analysis that we presented here concerning the CP estimation with correlation measures 
at finite temperature has also been considered in another study, where the entanglement of formation 
(EOF) and QD are the subjects of the investigation. As compared to the performance of EOF and 
QD on the same task, it is not possible to make an absolute comment on the success or failure of the 
measures considered in this work, due to its high dependence on the system parameters. However, 
relative comparisons can be made, for example in the case of nearest neighbors with 7 = 0.5 OMQC or 
in the case of next nearest neighbors with 7 = 0.001 WYSIM, outperforming both EOF and QD. 

We continue our discussion with the dependence of correlations on the distance between the spin pair 
for which they are calculated. It is known that bipartite entanglement is very short ranged in contrast 
to the QD. In Figure 4, we demonstrate our findings on the long-range behavior of WYSIM, MIN and 
OMQC for A = 0.75, 0.95, 1.05, 1.5, 7 = 0.001, 1 and kT = 0.1, 0.5. It can be seen from the plots 
that temperature have a diminishing effect on the correlations between distant spins. We also show that 
for 7 = 0.001, all measures decay to zero as the distance between spin pairs are increased. On the other 
hand, in the Ising limit of our model, 7 = 1, correlations between distant spin pairs settle at a finite value 
for A = 1.5, which is deep in the ferromagnetic (ordered) phase. In fact, in all cases considered here, 
long-range correlations persist longer in the ordered phase. 
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Figure 4. Long-range behavior of the thermal total and quantum eorrelations for 7 = 0.001 
and 7 = 1 at kT =0.1, 0.5. The eireles, squares, diamonds and triangles correspond to 
A = 0.75, A = 0.95, A = 1.05 and A = 1.5, respectively. 


0.6 


0.4 


0.21 


7=1 


7=1 

kT=0.1 


kT=0.1 _ 

* - A ■ - ■ A - A- ■ A - - * ■ • A - - * - • * - - ■ 



♦ 


♦. 




* < < : : a 1 ■ 1 ■ ■ 




10 0 


0 

10 


0.2 


A 



7=i 


7=1 

kT=0.5 


kT=0.5 

A 



♦ 


♦ *, 






A, 




'-A. 











._» t -- ml—tl—l 


0.3 


0.2 


0.1 


8 10 0 


8 10 


0.3 


0.2 \ 


0.1 


7 = 0.001 

kT=0.1 






7 = 0.001 

kT=0.1 


8 10 0 


1 < A 


0.3 


0.2 


0.1 


8 10 































Entropy 2014 , xx 


19 


4.2. Local Quantum Coherence 


As introduced in the eorrelations seetion, LQC quantifies the eoherenee eontained in a given density 
matrix with respeet to a fixed basis. The eoherenee measure that we adopt in this work is merely 
the WYSI. We start by the diseussion of single-spin eoherenee (eoherenee eontained in p when 
measuring ax), I{po, ctx) and its experimentally-friendly lower bound, I^{po, ax) in our system, whieh 
are presented in Figure 5. We plot these quantities for 7 = 0.5 and 7 = 1, whieh is the Ising model in 
transverse field. 

At this point, it is important to note that we have also ealeulated the eoherenee with respeet to different 
observables for both single- and two-spin eases; however, sinee they give no additional information about 
the system under eonsideration, we only show ax eoherenee as a representative. 

The left panel in Figure 5 shows the measures themselves, and the right panel displays their derivative. 
First, we reeognize the faet that the single-spin ax eoherenee and its lower bound deereases as the inverse 
field A inereases, meaning that in the ordered phase, a randomly ehosen spin eontains less eoherenee than 
it eontains the disordered phase. It ean be seen from Figure 5b,d that the derivatives of the measures are 
eapable of deteeting the presenee and the order of the CP by displaying a divergent behavior at Ac = 1. 
On the other hand, we see no non-trivial behavior at the faetorization field A/ ~ 1.1547 for 7 = 0.5, 
i.e., the eoherenee of a single spin behaves smoothly at this point. We do not expeet to see any effeet of 
faetorization for 7 = 1, sinee in this ease, the faetorizing field tends to infinity, A —)■ 00. It is remarkable 
and important that even the experimentally-friendly, simplified version of LQC, I^{po, ax), whieh does 
not require the full tomography of the state, ean spotlight the CP of the QPT. 


Figure 5. Single-spin -eoherenee for 7 = 0.5 (a) and 7 = 1 (c), along with its first 
derivative (with respeet to A) for 7 = 0.5 (b) and 7 = 1 (d), as a funetion of A. As the red 
solid line denotes the measure, the dashed blue line eorresponds to its simplified version. 
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Our results on two-spin LQC are presented in Figure 6, again for two different values of 7, 7 = 0.5 
and 7 = 1. Recall that in the two-spin case, we calculate the coherence contained in one of the 
subsystems locally. Therefore, the measurement operator only acts on the chosen subsystem, and the 
mathematical expression for this is given as / (pab, We have considered nearest neighbor spins 

in our discussion; however, similar results can also be obtained considering next-nearest neighbor spins. 
In the left panel, we can see the behavior of the measures and conclude that they seem very similar to that 
of the single-spin coherence. Both the original definition and its lower bound follow a decreasing trend 
with increasing inverse field. The derivatives of the measures, presented in the right panel, again signal 
the existence of the CP via a divergence in their derivatives, thus also giving the correct information about 
the order of the phase transition. The important difference from the single-spin case is the appearance 
of a small discontinuity in the original measure for 7 = 0.5 at A = 1.1547. As we have mentioned, this 
is the critical field in order to observe the factorized ground state of the system for the considered value 
of 7. While the WYSI can detect the occurrence of this phenomenon, the lower bound for WYSI does 
not get affected by it. It is rather striking that WYSI itself can signal the factorized ground state, since 
neither is it an entanglement measure for mixed states, nor have we considered the effects of symmetry 
breaking on the ground state of the system. 


Figure 6. Two-spin local cr^-coherence for 7 = 0.5 (a) and 7 = 1 (c), along with its first 
derivative (with respect to A) for 7 = 0.5 (b) and 7 = 1 (d), as a function of A. As the red 
solid line denotes the measure, the dashed blue line corresponds to its simplified version. 
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Now, we will have a deeper discussion of the occurrence of discontinuities in the derivatives of WYSI 
and WYSIM at A/. First of all, it is important to note that these discontinuities are present for all 
values of 7, not just for the 7 = 0.5 case. The explanation of such a behavior stems from the presence 
of the square root of the density matrix, y//^, in the definition of WYSI. The elements of are 
themselves discontinuous at the factorization point, resulting in a physical quantity depending on them 
to be also discontinuous at the same point. Only correlation measures other than WYSI that can identify 
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the factorization phenomenon in a single evaluation are entanglement measures, such as concurrence and 
the entanglement of formation, which is itself a function of concurrence. Interestingly, looking at the 
definition of concurrence, we see that it also involves and vanishes at A/ as a result of the fact that 
even the thermal ground state is factorized at this point. Furthermore, as soon as we drop the square root 
of the density matrix, for example as in the case of J^, we suddenly lose the information on the ground 
state factorization. We stress that the detection mechanism for the factorization point is not the same 
as the detection of the CP, since no thermodynamic quantity is discontinuous at this point. To be clear, 
we do not present an explicitly constructed correspondence between the discontinuities of the elements 
of and the emergence of the factorization point; instead, we point to a possible explanation of the 
results that we have obtained. 

In this last part of this section, we carry out the same investigation that we have done for the correlation 
measures in the previous section and look at the performance of LQC and its lower bound in correctly 
estimating the CP of the QPT and the factorization point. We have already explained that at temperatures 
low enough to see the effects of quantum fluctuations, it is possible to see the signatures of the CP, this 
time not through the divergences in the derivatives, but rather the extremum points near the CP. Similarly, 
at finite temperatures, the factorization point ceases to be a point and becomes a region of separability 
in the parameter space of the model. Finite temperature analysis is particularly important also for the 
reason that zero temperature is not attainable experimentally. 

In Figure 7a, we plot the estimated CP as pointed out by the LQC and its simplified version for 
7 = 0.5 as a function of temperature. We see that the experimentally-friendly lower bound remains 
quite accurate up to a temperature kT = 0.1. This is particularly important, since it proves that 
can be a strong candidate to detect the CP in experimental applications. Figure 7b shows the estimated 
factorization point calculated from the non-trivial behavior of a^, <Jy and coherence for 7 = 0.5 as 
a function of temperature. All three measures behave very similarly as the temperature increases. They 
stay extremely accurate until kT = 0.015 and start to deviate significantly from the actual point from 
kT = 0.02 

Figure 7. (a) The critical point estimated by single-spin cr2;-coherence (red line) and its 
simplified version (blue line) as a function of the temperature for 7 = 0.5. (b) The 
factorization field estimated by local two-spin cxa;-coherence (red line), cxj^-coherence (blue 
line) and cr^-coherence (green line) as a function of time for 7 = 0.5. 
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4.3. Local Quantum Uncertainty 

Lastly, we would like to diseuss the results that we obtain on LQU, whieh is nothing but the optimized 
version of the LQC over all possible observables. In Figure 8, we plot LQU and its derivative for 
7 = 0.5 and 7 = 1 as a funetion of A. The measure itself has a small asymmetry around the peak 
point near A = 1, settling to a higher value in the ordered phase. Through the divergenee at Ac = 1 
and finite diseontinuity at A/ = 1.1547, in its derivative, we ean spot the presenee of the CP and the 
faetorization point, respeetively. However, in Figure 8a,e, there are two extremum points, resulting in 
diseontinuities in the derivatives, at the Hamiltonian parameter values where the system exhibits no 
non-trivial behavior. The reason behind these non-analytieities in the derivative of LQU is, in faet, due 
to the optimization proeedure involved in the ealeulation of LQU. An abrupt ehange in the optimizing 
observable eauses sueh extremum points in the measure itself. Speeifieally, in this ease, two extremum 
points of LQU exaetly eorrespond to the values of A, where the optimizing observable ehanges from 
to ax- Correlation measures involving optimization proeedures might sometimes eause sueh ambiguities 
and display finite diseontinuities that are not rooted in the elements of the redueed density matrix. 
Therefore, when dealing with sueh measures, one should be eareful not to eonfuse a non-analytieity 
resulting from the optimization with a QPT. 

Figure 8. Two-spin loeal quantum uneertainty for 7 = 0.5 (a) and 7 = 1 (c), along with its 

first derivative (with respeet to A) for 7 = 0.5 (b) and 7 = 1 (d), as a funetion of A. 




5. Conclusions 

In this work, we have presented a review of the behavior of eorrelations in quantum spin ehains 
and their relation to the phenomena of faetorization and QPTs. We have provided an overview of the 
signifieant papers in the field, summarizing their main results. We have also ineluded an extensive survey 
of our eontributions to related problems. 
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In particular, we have firstly studied the thermal quantum and total eorrelations in the XY ehain. 
There are four different measures eonsidered in our investigation, introdueed in Seetion 3, namely 
MIN, WYSIM, OMQC and eoneurrenee. We have analyzed the behavior of these eorrelations near 
the CP of the QPT and seen that all four measures are able to deteet the presenee of the QPT at zero 
temperature. As we inerease the temperature, measures start to lose their non-trivial behavior near the 
CP Therefore, divergent derivatives are replaeed with extrema points. However, looking at the loeation 
of these extrema, it is possible to determine an estimated CP at finite temperatures. In faet, we have 
used this method to eompare the aeeuraey of the estimation performanee of the eonsidered measures at 
finite temperatures. We have seen that the performanee of the measures in eorreetly estimating the CP 
of the QPTs is strongly dependent on the Hamiltonian parameters. For example, while OMQC has good 
performanee in estimating the CP at finite temperature for 7 = 0.5, its aeeuraey signifieantly deereases 
for 7 = 0.001. Similar results have also been obtained for seeond nearest neighbors, but sinee they do 
not provide any new information about the system, we have not ineluded them in our review. We have 
also explored the long-range behavior of the eorrelations and seen that in the ordered (ferromagnetie) 
phase, they remain non-zero for very large spin separations. 

Seeond, we analyze the behavior of the eoherenee measure WYSI in the XY ehain. We have seen that 
even if we look at the eoherenee for a single spin, we ean identify the QPT by the divergenee in the 
first derivative. Furthermore, the experimentally-friendly simplified version of the eoherenee measure 
ean also deteet the QPT, whieh is an important result for possible experimental applieations. On the 
other hand, both of them fail to reeognize the faetorized ground state. Moving on to the two-spin 
loeal eoherenee, we again observe that the loeation of the CP ean be spotlighted by the divergenee in the 
measure. However, even though the simplified measure again does not feel the presenee of the FP, LQC 
is able to pinpoint the position of the FP. Therefore, we have eoneluded that the diseontinuous behavior 
of a eorrelation measure at the FP stems from the non-trivial behavior of the elements of the square root 
of the density matrix. It is important to note that using different observables from does not provide any 
additional insight about the system under eonsideration. We have then turned our attention to the LQC, 
whieh is the optimized version of the LQC over all sets of observables. The divergenee of the derivative 
of LQU is again present at the CP pinpointing the QPT. However, we see sharp extrema points, whieh do 
not eorrespond to any kind of non-triviality of the XY model. A detailed analysis of these points revealed 
that the extrema points are in faet eaused by the ehange in the optimizing observable in the definition of 
LQC. Thus, it is important to be eareful about the origin of non-analytieities showing-up in the behavior 
of quantum eorrelations before reaehing eonelusions about eritieality and faetorization. 
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